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ABSTRACT 

Halos formed in the standard A cold dark matter framework should follow an universal 
mass density profile and fit a well defined mass-concentration relation. Lensing analyses 
of clusters with a large Einstein radius seem to contradict this scenario, with the massive 
cluster Abell 1689 being often claimed as a notable example of a highly over-concentrated 
halo. Shape and orientation biases in lensing studies might be at the basis of this disagree- 
ment between theory and observations. We developed a method for a full three-dimensional 
analysis of strong and weak lensing data. Surface density maps estimated from lensing are 
de-projected to infer the actual triaxial structure of the cluster, whose mass distribution is 
approximated as an ellipsoidal Navarro-Frenk- White halo with arbitrary orientation. Inver- 
sion is performed under competing a priori assumptions, integrated in the method thanks 
to Bayesian statistics. We applied the method to A1689. Whatever the considered priors on 
shape and orientation, both weak and strong lensing analyses found the halo to be slightly 
over-concentrated but still consistent with theoretical predictions. We found some evidence 
for a mildly triaxial lens (minor to major axis ratio ^ 0.5 ± 0.2) with the major axis ori- 
entated along the line of sight. Exploiting priors from A^-body simulations, we found mass 
-^^200 = (1.3 ±0.4) X IO^^A/q and concentration C200 = 10 ±3 for the weak lensing analysis 
of Subaru data, A/20o = (1-7 ± 0.3) x 10^^ Af© and C200 = 6.1 ± 0.9 for the strong lensing 
analysis of multiple image systems, and A/200 — (1-3 ± 0.2) x IO^^Mq and C200 = 7.3 ± 0.8 
for the combined weak plus strong analysis. 

Key words: galaxies: clusters: general - cosmology: observations - gravitational lensing - 
methods: statistical - galaxies: clusters: individual: Abell 1689 



1 INTRODUCTION 

Clusters of galaxies are the most recent bound structures to 
form in a hierarchical cold dark matter model with a cosmo- 
logical constant (ACDM). A precise determination of their mass 
and shape offer important clues to the understanding of the as- 
sembly process of structure in the universe. A-body simula- 
tions are successful in making detailed theoretical predic tions on 
dark matter halo properties dNavarro et al. 1997; Bull ock et al.l 
I2OOII : lDiemandetij]|2004l : iDuffv et alj 12008. : Prada et aTTEoilh . 
and lensing measurements are supposed to provide precise and 
model independent mass measurements. W hereas the universal 
Nava rro-Freank-White (NEW) density profile dNavarro et al.l 199^ . 

reproduces many characteristics of massive lenses, some 
puzzling disagreement between the ACDM theoretical frame- 
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work and measurements still persists. One possible conflict is 
the detectio n of extremely large Einstein radii in massive lens- 
ing clusters jSroadhurst & BarkanalbOOSHSadeh & Rephaelil2008l : 
lOguri & Blandfordll2009l: IZitrin et al.ll2011allbl) . Haloes should be 



over-concentrated to fit the data. 

Concentrations of massive galaxy clusters are a crucial probe 
of the mean density of the universe at relatively late epochs. The 
concentration parameter measures the halo central density and 
should be related to its v irial mass, with the c oncentration decreas- 
ing gradually with mass dBullock et al.l200ll) . However, cluster ob- 
servations have yet to firmly confirm this correlation. 

The observed concentration-mass relation for galaxy clus- 
ters has a slope consistent with what found in simulations, 
thoug h the normalization factor is higher dComerford & NataraianI 
l2007h . Lensing concentrati ons appear to be systematical ly larger 
than X-ray concentrations dComerford & Nataraianll2007k . A sim- 
ilar, though less pronounced, effect is also fo und in simulations 
dHennawi et alj l2007l : iMeneghetti et alj l201(t) . which show that 
massive lensing clusters are usually elongated along the line of 
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sight. lOguri & BlandfordI ( l2009h showed that the larger the Ein- 
stein radius, the larger the over-concentration problem, with clus- 
ters looking more massive and concentrated caused by the orienta- 
tion bias. 

The observational picture from lensing analyses is still not 
clear. iBroadhurst et al.l (jioOS) derived lens distortion and magni- 
fication of four nearly relaxed high-mass clusters, inferring sig- 
nificantly high concentrations. The analysis of som e clusters with 
very large Einstein radii (Zitrin et al. 20091 20 11 ah als o showed a 
strong over-concentration problem. Ogu ri et al. [ (l2009b found that 
the data from a sample of ten clusters with strong and weak lensing 
features were highly inconsistent with the predicted concentration 
parameters, even including a 50 per cent enhancement to account 
for the le nsing bias. Strong lensing features in the MACS cluster 
sample (Zitrin et^l2011a : lMeneghetti et al.ll201lh confirmed the 
over-concentration problem. 

On the other ha nd, out of 30 X-ray clusters with significant 
weak lensing signal, lOkabe et al.l feOlOl) found that the 19 clus- 
ters that were well fitted by a NEW profile showed a correlation in 
the c{M) relation which is marginally compatible with predictions 
for both slope and normalization. Einally, weak lensing analyses 
of stacked clust ers of lesser mass do not exhibit the high conc en- 
tration problem (Johnston et al. 2007; Mandelbaum et al.'2008|), in 
agreement with the findings of Oguri & Blandford (2009). 

Several effects can play a role. Strong lensing clusters tend 
to prefer entially sample the high-mass end of the cluster mass 
function jComerford & NataraianI l2007h . Lensing concentrations 
can be also inflated due to substructures close to the line of 
sight. Eurthermore, without adequate colour information, the weak- 
lensing signal can be diluted particularly toward the cluster cen- 
tre by the inclusion of unlensed cluster members in the definition 
of background galaxy samples (Broadhurst et al. 2005), leading 
to an underestimation of the concentration (iLimousin et al.ll2007l : 
lUmetsu & Broadhurstll2008h and of the central cluster mass exam- 
ined at high overdensities ( Okabe et^l2010h 

Disagreement between theory and observation might be ex- 
plained by some orientation and shape biases. Triaxial halos can 
be m uch more efficient lenses than their more spherical counter- 
parts dOguri & BlandfordI 120091) with the strongest lenses in the 
imiverse expected to be a highly biased population preferentially 
orientated along the line of sight and with high levels of triaxiality. 
On the other hand, Hamana et al. ( 2009) carried out a spectroscopic 
follow-up survey of weak- lensing selected cluster candidates of 
iMivazaki et"ZI I2002! ). and found no evidence of lensing selection 
bias in mass estimates compared to X-ray and optical techniques. 

Only a few works have tried to infer mass and concen- 
tration from lensin g data without assuming spherical symmetry. 
lOguri et^ ( l2005h studied the mass profile of A1689 exploiting 
both strong and weak lensing measurements. They showed that the 
intrinsic model parameters are weakly constrained (mainly because 
of the lack of information along the line of sight), and that A 1689 
lensing observations might be compatible with predictions assum- 
ing a major axis of the hal o aligned with the lin e of sight. Triaxiality 
issues were addressed bv lCorless et al.l ( l2009h , who put weak lens- 
ing constraints on three strong lensing clusters exploiting Bayesian 
statistical methods. They found that the large errors that accom- 
pany triaxial parameter estimates can make observations compat- 
ible, even if marginally, with theoretical predictions. In fact, ne- 
glecting halo triaxiality can lead to over- and under-estimates of up 
to 50 per c ent and a factor of 2 in halo mass and conc entration, re- 
spectively dCorless & King|2007l ). ISereno et alj ( l2010l) investigated 
a sample of 10 strong lensing clusters considering which intrin- 



sic shape and orientation the lensing halos should have to account 
for both theoretical predictions and observations. They found that 
nearly one half of the clusters seemed to be composed of outliers of 
the mass-concentration relation, whereas the second half supported 
expectations of A'^-body simulations which prefer mildly triaxial 
lensing clusters with a strong orientation bias. The strong lensing 
analysis of AC 1 14 also supported that mapping effects can play a 
role in the estimate of the concentration (iSereno et al.l2O10l) . 

Here we apply a novel lensing inversion method to infer in- 
trinsic cluster parameters. We exploit both weak and strong lens- 
ing data through a full two-dimensional (2D) analysis and then de- 
project the lensing derived surface mass density with a Bayesian 
approach. The method is applied to Abell 1689 (A1689), a very 
massive cluster at redshift z = 0.183 with one of the largest Ein- 
stein radius observed to date, ^ 45" for a s ource redshift Zs = 1 
dBroadhurst et al.|[2005ULimousin et al.ll2007i) . Despite being thor- 
oughly investigated it is still unclear if and how much this cluster 
is over-concentrated. Independent strong lensing analyses of the 
inner regions provided a co nsistent picture of the mass distribu- 
tion in the inner < 3()0 kpc i Broadhurst et al.ll2005l : lHalkola et"ZI 
I2OO6I : iLimousin etZll2007l : ICoe et alj|2010h , but results are not 
fully in agreement on the larger virial scale. Different weak 
lensing an alyses of either Subaru (Umetsu & Broa dhurst! l2008l : 
lUmetsuet al. 2009). ESO/MPG (Coriess et al. 2009), or CEH12K 
I Limousin et al.ll2007l) suggest somewhat different degrees of con- 
centration, C200- Subaru and CEHT results are marginally consis- 
tent: C 200 ~ 10.1 ± 0.8(=b2.2 s ystematic) in Umetsu & Broadhurst! 
( !2008!) : C200 = 7.6 ± 1.6 in [Limousin et al.! (!2007^. The broad 
conclusion is that weak lensing investiga tions point out t o a more 
concentrated cluster than strong lensing (!Coe et al.!!2OI0l) . but the 
quantitative assessment of such tension is still to be firmly ascer- 
tained. 

One way to possibly solve the overconcentration problem in 
A 1 689 is to prope rly consider the triaxial structure of the cluster. 
!Oguri et al.! J2005!) first noted that measured concentration in lens- 
ing analyses of A1689 would be marginally in agreement with the- 
oretical predictio n if the maj or axis of the cluster is aligned with 
the line of sight. iDe Eilippis et al.l ( |2005|) and ISereno et al.! ( I2OO6I) 
combined X-ray and Sunyaev-Zel'dovich (SZ) observations to in- 
fer the elongation of the hot intra-cluster medium without assuming 
hydrostatic equilibrium and found that A 1689 was quite spherical 
within the errors. A triaxial weak lensing analysis performed in 
ICorless et al.! ( l2009l) revealed that realistic intrinsic parameter esti- 
mates in A 1689 must be accompanied by larger errors due to our 
limited knowledge of the halo shape and orientation. 

The h ypothesis of h ydrostatic equilibrium in A 1689 was ex- 
ploited by iLemze et al.l ( I2009) . who combined Chandra X-ray 
brightness measurements and joint strong/weak lensing measure- 
ments of the projected cluster mass profile by IBroadhurst et al.l 
( I2OO5I) to derive a temperature profile that is consistent with hy- 
drostatic equilibrium. It was found that the resulting equilibrium 
temp erature exceeds the observed temperature by 30 per cent at all 
radii . iLemze et al. I( l200l suggested that the existence of significant 
cool and dense gas components is the source of this temperature 
discrepancy (Kawahara et al. 2007). 

On the basis of a d etailed X-ray analysis of deep Chandra 
observations, IPeng et al.l Bm% showed that the hydrostatic mass 
within the central 200 kpc /i~^region of A1689 is about 30-50 per 
cent lower than lensing-based mass estimates obtained under the 
assumption of spherical symmetry. They also showed that it is un- 
likely that dense and cool gas clumps alone can cause such a strong 
bias in the X-ray temperature determinations, since the space fA\- 
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ing factor for such components should be 70-90 per cent. They 
suggested instead that a prolate ellipsoidal configuration, aligned 
with the line of sight, could solve the central mass discrepancy as- 
suming an axis ratio of 0.6 for the gas distribution, which however 
overesti mates the pr ojected total mass at large radii by ~ 40 per 
cent. Mo randi et alj (2011) confirmed that if the halo of A 1689 is 
very elongated along the line of sight, strong lensing and X-ray ob- 
servations are compatible with theoretical predictions on the c(M) 
relation. 

The level of hydrostatic equilibrium in A1689 has 

been recent ly explored by .Kaw aharad a et al.l and 

Molnar et alj koipi) . Dee p mosaic Suzaku X-ray observations by 
Kawaharada et alj ilOld) . covering the entire cluster, revealed 
anisotropic temperature and entropy distributions in cluster out- 
skirts correlated with large-scale structure of galaxies. Three out 
of four quadrants of the mosaic image, in contact with low density 
void environments, have low gas te mperatures and entropies , 
deviating from hydrostatic equilibrium. ^ Kawaharada et al.l bOlOl) 
showed that the thermal gas pressure within < 1.3 Mpc (about half 
the virial radius) is at most about half of the equilibrium pressure 
to balance the gravity predicted by gravitational lensing under the 
hypothesis of spherical symmetry, and ~ 30 — 40 per cent around 
the virial radius. They suggested that additional sources of pressure 
support, such as bulk and/or turbulent motions, are required in the 
cent ral and outskir t regions. 

iMolnar et al. challenged the assumption of hydrostatic 

equilibrium in the core region of A1689. By analysing a sample of 
high-mass clusters with M = (1 — 2) x 10^^ AIq drawn from high- 
resolution cosmological simulations, they found in the core region 
of relaxed clusters, selected based on their smooth spherically av- 
eraged density profile and no sign of recent major merger events, 
significant non-thermal pressure support of ~ 20 per cent from 
subs onic random gas motions within 1-10 per cent of the virial ra- 
dius jYounger & Brvanl2007l:lLau et al.l2009i) . They suggested that 
discrepancies might be due in part to non-thermal pressure support. 

Here, we reconsider the estimate of the concentration in 
A 1689 from lensing analyses renouncing the spherical hypothesis. 
The paper is organised as follows. In Sec.|2] we review the prop- 
erties of the three-dimensional NFW profile and how it projects 
on the plane of the sky. Theoretical predictions derived in TV-body 
simulations are discussed in Sec. [3] Fitting results to strong and 
weak lensing data can be found in Sec.|4l The lensing inversion is 
discussed in Sec.|5] whereas Sec.|6]focuses on results. Section|7]is 
devoted to the conclusions. Throughout the paper, we assume a flat 
ACDM cosmology with density parameters JIm ~ 0.3, JIa = 0.7 
and Hubble constant Ho = lOO/i km s"^Mpc"\ h = 0.7. At the 
A1689 distance, 1" corresponds to 3.08 kpc. 



2 BASICS ON TRIAXIAL NFW HALO 

High resolution A^-body simulations have shown that the density 
distribution of d ark matter halos is successf ully described as a NFW 
density profile dNavarro et al.lll996l Il997h . whose 3D distribution 
follows 



Pnfw ~ 



{r/rs){l + r/rsy 



(1) 



where ps is the characteristic density and Vs is the characteristic 
length scale. Halos are expected to be aspherical and their profiles 
can be accura tely described by concentric triaxial ellipsoids with 
aligned axes tfing & Sutdl2002h . The spherical radius r has to be 



replaced with an ellipsoidal radial variable. The minor (intermedi- 
ate) to major axial ratio is denoted as gi (g2)with0 < gi ^ §2 ^ 1; 
we also use the inverse ratios, < = 1/qi ^ 1. 

Different definitions of parameters for spherically averaged 
profiles can play a role when comparing observations to predictions 

iBroadhurst & Barkana 2008). Here, we follow Corless & Kiri^ 
20071) . who generalised NFW parameter definitions to ellipsoidal 
halos. The radius r2oo is such that the mean density contained 
within an ellipsoid of semi-major axis r2oo is A = 200 times the 
critical density at the halo redshift, pcr', the corresponding concen- 
tration is C200 = 7'2oo/7's. M200 is the mass within the ellipsoid 
of semi-major axis r2oo, M200 = (8007r/3)gi(j2?'2oo/5cr. Such de- 
fined C200 and M200 have small deviations with respect to the pa- 
rameters computed fitting spherically averaged density profiles, as 
done in A^'-body simulations. The only caveat is that the spheri- 
cal mass obtained in simulations is si gnificantly less than th e ellip- 
soidal M200 for extreme axial ratios ( Corless & King|2007t) . How- 
ever, since the dependence of the concentration on the mass is quite 
weak, see Eq. il3l . effects on our analysis are negligible. 

Three Euler's angles, 0, ip and relate the intrinsic to the 
observer's coordinate system; the angle 9 quantifies the inclination 
of the major axis with respect to the line of sight. When viewed 
from an arbitrary direction, quantities constant on similar ellipsoids 
project themselves on similar ellipses (iStark.l977i) . 



2.1 Projected density 

In gravitational lensing studies, the projected map of interest is 
the surface mass density, expressed in terms of the convergence 
K, i.e., in units of the critical surface mass density for lensing, 
Scr = {c^ Db) / {4tvG Dd Dds), where Ds, -Dd and D^s are the 
source, the lens and the lens-source angular diameter distances re- 
spectively. The projected surface density has the same functional 
form as a spherically symmetric halo. The axial ratio of the major 
to the minor a xis of the obser ved projected isophotes, ep 1), can 
be written as ( lBinggelilll980l) . 



ep 



where j, k and / are defined as 



(2) 



2 2 • 2 

J = 6162 sm I 



2 2 2 2 2 2 

61 cos 6 COS + 62 COS 6 sin ip, (3) 



k = (ci — 62) sin (/p cos cos ^, 
I — ei sin^ ((9 + 62 cos'^ (p. 



(4) 
(5) 

We also use the ellipticity e = 1 — 1/ep. 

The projected NFW density profile can be described in terms 
of the strength of the lens ks and of the projected length scale Vsp, 
i.e., the two param eters directly inferred by fitting lensing maps 
al.l2OI0h . The observed scale length rsp is t he projection 
on the plane of the sky of the cluster intrinsic length ( IStark|[l977l : 
ISereno 2007), 

f £^ 



ereno et 



The observed scale length VsP is t he projection 



where / is a function of the cluster shape and orientation, 

2 2 2 2 2 2 2 

f = ei sin 9 sin <^ + 62 sin 9 cos (p + cos 0; 



(6) 



(7) 



the parameter eA quantifies the elongation of the triaxial ellipsoid 
along the line of sight. 



\ 1/2 
, ep \ 3/4 

,6162 / 



(8) 
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ga represents the ratio between the major axis of the projected el- 
lipse in the plane of the sky and the size of the ellipsoid along the 
line of sight. When eA < 1, the cluster is more elongated along the 
line of sight than wide in the plane of the sky, i.e., the smaller the 
fiA parameter, the larger the elongation along the line of sight. 

The relation between the projected par ameters, fcs and r^p , 
and the intrinsic parameters can be written as jSereno et al.llOlOh . 



E7 ./gCO 
cr X Kg — — ^^Ps^sPi 



where as usual ps — Scpcr, 
200 C200 



3 ln(l + C20o) - C20o/(l + C20o) ' 

and /gco = v^ep/eA- The mass M200 can be expressed as 



M200 = ^ X 200pcr X (c2oorsp)^%^. 
3 ep' 



3 THEORETICAL PREDICTIONS 



(9) 



(10) 



(11) 



iV-body simul a tions jNeto et al.l |20 07|: iMaccio et alj I2OO8I : 
iGao et all |2008| ; iDuffv et al.l I2OO8: Pra daet alj|201ll) have pro- 
vided a quite detailed picture of the expected properties of dark 
matter halos. Results may depend on the overall normalization 
of the power spectrum, t he mass resolution and the simulation 
volume ( Prada et al.ll201 ll) . The dependence of halo concentration 
on mass and redshift can be adequately described by a power law, 

c = A{M/Mpi.otf{l + zf . (12) 

Here, we follow iDuffv et alj bOOSh , who used the cosmologi- 
cal parameters from WMAP5 and found {A, B, C} = {5.71 ± 
0.12, -.084 ± 0.006, -0.47 ± 0.04} for a pivotal mass Mpivot = 
2 X 10^^Mq//i in the redshift range — 2 for their full sample of 
clusters. The scatter in the concentration about the median c(M ) 
relation has been shown to be lognormal. 



p{lnc\M) 



1 



exp 



Inc - In c(M) 



(13) 



with a dispersion o ^(lo^^ n C200) ~ 0.15 for a full sample of clusters 
jPuffv et al J2008t). 

|jing&Sutd ( l2002l) investigated the probability distribution of 



intrinsic axial ratios and proposed an universal approximating for- 
mula for the distribution of minor to major axis ratios. 



p((ji) oc exp 



(gi -<lt^/rqif 



where — 0.54, a. 



2a2 

= 0.113 and 

0.7 



(14) 



(15) 



with the characteristic nonlinear mass at redshift z and Afvir 
the virial mass. The conditional probability for §2 goes as 



p{qi/q2\qi) 



2(1 



2gi/i32 



(16) 



for qi/q2 ^ ^'min = max[gi,0.5], whereas is null otherwise. The 
lensing population has n early the same triaxia lily distribution as the 
total cluster population (" Hermawi et al.ll2007h . 

For comparison we also consider a flat distribution for the ax- 
ial ratio in the range qtmin < qi ^ 1 and ^ 52 5; 1. Probabilities 
are defined such that the marginalized probability p(qi ) — const. 



and the conditional probability p{q2\qi) ~ const. The probabili- 
ties can then be expressed as 



P(9l) = 1/(1 - "Jmin) 

for Qmin < gi ^ 1 and zero otherwise; 

p{q2\qi) = (1 - qiy^ 



(17) 



(18) 



for q2 J? qi and zero otherwise. The resulting probability for §2 is 
then 



P(92) 



1 - gn 



In 



1 — gmiii 
1-52 



(19) 



A flat distribution allows also for very triaxial clusters (gi < §2 ^ 
1), which are preferentially excluded by A'^-body simulations. 

Finall y, semi-analytical lOguri & BlandfordI l2009b and nu- 
merical ( iHennawi et alj|2007l) investigations showed a large ten- 
dency for lensing clusters to b e aligned with the li ne of sight. Such 
condition can be expressed as jCorless et al.ll2009l) 



p(cos 6) oc exp 



(cose - 1)^ 



2ai 



(20) 



with ae = 0.115. For comparison, we also consider a population 
of randomly oriented clusters with 

p{cose) = l (21) 
forO ^ cose s; 1. 

4 LENSING ANALYSIS 

In this section, we fit strong and weak lensing data with a projected 
ellipsoidal NFW halo. 

4.1 Weak lensing 

A1689 has been object of very detailed weak lensing (WL) studies. 
Here, we rely on previo us work by Umetsu & Broadhurst (2008|) 
and Umets u et al.l ^20091) on wide-field Subaru data. They obtained 
the two-dimensional convergence map reconstructed from a reg- 
ularised maximum-likelihood combination of the joint measure- 
ments of the weak-lensing distortion and magnification bias effects 
on background red galaxies, the combination of which can break 
the mass-sheet degeneracy inherent i n all reconstruction m ethods 
based on distortion information alone l lBroadhurst et al Tl995il . The 
map is given on 21 x 17(= 357) grid pixels covering a field of 
^ 30' X 24'. Such a map probes the cluster halo up to and be- 
yond the virial radius (~ 16') with the accuracy of the pixel size 
(— 1.4'). The map we considered was obtained with a pure WL 
(shear plu s magnification) recons t ructio n method, described as "2D 
MEM" in lUmetsu & BroadhursJ ( l2008h . 

A full two-dimensional weak lensing anal ysis is cru cial for 
comparison with theoretical predictions (Oguri et al ] |2010l) . Given 
the measured convergen ce map fcpbs, the w eak lensing Xwl func- 
tion can be expressed as dOguri et al.ll2005h . 



2 

Xwl 



= [fcobs(rO - fc(rO] (F-i)^ . [fcobs(r,) - fc(r,)] (22) 



where is the inverse of the pixel-pixel covariance matrix. The 
corresponding likelihood is £wl oc exp(— Xwl/2)- The x^(— 
319) for the best fitting model (found with a downhill simplex al- 

i ;orith m) is very simila r to w hat found in lUmetsu & BroadhurstI 
20081) or lUmetsuet"ZI ( l2009h . The slight improvement is due to 
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Component 


mass scale 




62,0 


e 




length scale 






y ) 


\ ) 








NFW 


0.35(0.37) ±0.04 
H") 


0.0(0.4) ±0.6 


Strong Lensing 
-4.6(-4.4)±0.6 


0.29(0.27) ±0.04 


-26(-24) ±5 


r-s-p 

400(370) ± 70 


p-Jaffe 


15(13) ±5 


-80(-70) ±30 


80(70) ±40 


0.88(0.88) ±0.06 


-142(-143) ±2 


43(38) ±23 


NFW 


0.62(0.53) ±0.24 


-6(-6) ±7 


Weak Lensing 
-5(-7)±5 


0.17(0.21) ±0.12 


30(30) ±30 


210(230) ±70 




ks 

0.41(0.42) ±0.03 




Weak and Strong Lensing 






NFW 


-0.1(0.1) ±0.6 


-4.6(-4.5)±0.6 


0.30(0.29) ±0.04 


-22(-21) ±4 


310(300) ±30 



Table 1. Projected model parameters inferred with the lensing analysis. The orientation angle df is measured north over east and r^p (r^p) is the projected 
length scale (core radius) for the NFW (p-Jaffe) profile. Convergence is normalised to a reference source redshift (2a = 2.0). The truncation radius of the 
p-Jaffe component was kept fixed to 500 k pc. For the W L plus SL analysis, we report only the results for the main halo. Central location and dispersion of the 
PDF are computed as biweight estimators teeers et al.ll990h . Maximum likelihood values are reported in parentheses. 
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Figure 1. Marginalised posterior probability density functions for the projected parameters of the NFW halo. Panels from the left to the right are for the biimed 
PDF of ks, Tgp, and e, respectively. Shadowed and white bins are for the weak and strong lensing results, respectively. 



the use of an elliptical model. We run four Markov chains of 50000 
samples each after an early parameter diagonalization. Chain con- 
vergence was reached. For each series of chains run in the present 
paper, we checked that the standard var(chain mean)/mean(chain 
var) indicator was less than 1.2. Results are summarised in Table[T] 

4.2 Strong lensing analysis 

The HST observations of Abell 1689 are among the deepest to 
date for any galaxy cluster. The cluster has been observed in 
four filters ((7475 , rg25 , 2775, and zgso passbands) with ACS 
during 20 orbits of HS T and eight more bands from the ground 
teroadhurst et al.l2005h . The core of A1689 is very rich in detected 
lensi ng features: 135 candidate images of 42 galaxies were identi- 
fied jLimousin et al.ll2007l ; ICoe et al" We considered only 
systems with confirmed spectroscopic redshifts, which left us with 
80 confirmed images in 24 systemo 

We employed a parametric lensing analysi s method based on 
the gravlens kernel evaluation of xIl values jKeetodliooil lbh. 
As usual, strong lensing (SL) likelihood values were computed by 
comparing observed to predicted image positions. 

Since we are interested in constraining the total mass distribu- 
tion, we tried to keep the model as simple as possible. For an ideal 

1 Systems 1, 2, 4, 5, 6, 7, 10, 11, 12, 14, 15, 17, 18, 19, 22, 24, 29, 30, 31, 
32, 33, 35, 36, 40 according to the notation in lLimousin et alj )2007l) . 



comparison with theoretical predictions, we should model the total 
matter distribution in the cluster as a single NFW halo. However, 
parametric strong lensing models for A1689 need at least a second 
matter clump at northwest together with the main halo to be effec- 
tive. We did not try to model separately the halo of the brightest 
cluster galaxy (BCG) for two main reasons. Firstly, we were inter- 
ested in total matter distribution and we did not try to distinguish 
between baryons and dar k matter. This requi re the modelling of the 
intra-cluster medium too jSereno et al.l2010l) . Secondly, in order to 
effectively differentiate the BCG contribution from the main dark 
matter halo , we n eed dynamical data of the inner velocity profile 
(Sandetali2008l) . 

We could parametrize the inner cluster profile as a main 
NFW halo plus an additional north-west component, modelled as 
a pseudo-Jaffe profile (Keeton 2001b). As for weak lensing, we 
explored the parameter space by running four Markov chains of 
50000 samples each. Parameters were diagonalized after an early 
run. Results are summarised in Table[T] The average resi dual in im- 
age fitting was < 2.5" in line with previous modelings dCoe et al] 
[2010, see table 1). 

Strong lensing parametric models may be oversimplified. A 
fitting model with the accuracy of the HST astrometric resolution 
should resolve substructure a pproximately 25 kpc across within the 
central ~ 400 kpc diameter dCoe et al]|2010 '). Simple models can- 
not account for such mass complexities. A coarser positional error 
allowed us to perform the lensing analysis without adding too many 
parameters and, at the same time, avoiding the explored region in 
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0.0 0.5 1.0 1.5 2.0 

Figure 2. Contour plot of the marginalised PDF for the strength fcs and the 
projected radius r^p. Contours are plotted at fraction values exp(— 2.3/2), 
exp(— 6.17/2), and cxp(— 11.8/2) of the maximum, which denote confi- 
dence limit region of 1, 2 and Scr in a maximum likelihood framework, re- 
spectively. The shadowed regions are for the WL results, the thick contours 
for the SL. The inset zooms in for the SL region. 



parameter space being overly confined dSand et al.l2008h . We took 
the following approach. First, we found a global minimum of the 
likelihood function (xsl ~ 3231 for a positional uncertainty of 
/S.9 — 0.4"). Since all images were provided with the same er- 
ror, the chosen uncertainty does not affect the position of the mini- 
mum. Then, we rescaled the positional eiTor such that the reduced 
Xsl ~ 1- Such rescaled error was then used in the Markov chain 
analysis. 

Comparing weak to strong lensing results can be tricky. Strong 
lensing focuses on the cluster core and extrapolation well beyond 
the observational radius covered by multiple images is needed to 
infer global properties of the halo. On the other side, weak lensing 
observations provide a good coverage of the virial radius but are 
not accurate enough in the central regions to infer the inner slope 
of the density profile. Assuming the same parametric model for 
both strong and weak lensing assures us of interpreting the data co- 
herently in the same framework, but we have still to check for con- 
sistency. Projected parameters inferred from either weak or strong 
lensing analysis of A 1689 turn out to be in good agreement, see 
Table[T]and Figs.[T][2l Marginalised posterior distributions for the 
halo strength, projected radius and ellipticity are plotted in Figs.[T] 
Even if the strong lensing analysis prefers smaller (larger) values 
for fcs (»'sP, respectively), distributions are still compatible, as also 
enlightened by comparing the marginalised distribution p(fcs, Vsp), 
see Fig.|2] The SL constraint on ks is much tighter than the WL one, 
whereas dispersions in the probability density functions (PDFs) for 
the projected radius are similar. On the other hand, the weak lens- 
ing constraint on ellipticity is not very strong, since only very large 
values of e are ruled out. Being the WL PDF for ellipticity quite 
flat, the good agreement with the SL result is not very significant. 



5 LENSING INVERSION 

Deprojecting a surface density map to infer the intrinsic 3D shape 
is an under-constrained astronomical problem. Let us consider the 
ellipsoidal NFW halo. The 3D halo feature six parameters: the den- 
sity profile is described by two parameters (M200 and C200), the 
shape by two axial ratios (qi and g2), the orientation by two an- 
gles, 9 and (fl On the other hand, from lensing observations, we 
only get three constraints relating the intrinsic parameters: the lens- 
ing strength ks, the projected radius Vsp, and the ellipticity e. Even 
combining multi-wavelength observations, from X-ray to optical 
to radio band, one can only constrain the elongation of the cluster 
along the line of sight ( Sereno 2007,) . 

One wa y to address such shortcomings is to ex ploit some prior 
information joguri et al.l2005l : ICorless et al.|2069l) . The Bayes the- 
orem states that 

p(P|d) cx £(P|d)p(P), (23) 

where p(P|d) is the posterior probability of the parameters P 
given the data d, £(P\d) is the likelihood of the data given the 
model parameters and p(P) is the prior probability distribution for 
the model p arameters. The stat istical appro ach de s cribed above is 
the same of lOguri et al] ( l2005h and .Corless et al.l ( l2009h . but our 
method differs in some i mportant points as fa r as the strong lens- 
ing analysis is concerned. Corless et al. only analysed weak 
lensing data, whereas Oguri et al. (2005) used a circularly aver- 
aged profile of the convergence map of the inner regions to put 
strong lensing constraints. In the present paper, we have performed 
a parametric analysis of the multiple image systems so that we can 
employ a strong lensing likelihood that is fully triaxial. 

The likelihood for the intrinsic parameters is the same as in 
Sec. 14.21 and 14. 1 I for strong or weak lensing analyses, respectively. 
In fact, the three projected parameters can be expressed as functions 
of the intrinsic ones, see Sec.|2] 

We considered several priors, see Sec.|3] For the axial ratios qi 
and q2, we considered either the TV-body predictions in Eqs. ( 1141 - 
I16t or a flat distribution, see Eqs. ( I17II18) . We always put a lower 
bound qi ^ 0.05. For the alignment angle 6, we considered either 
the biased distribution for p(6) in Eq. J20b or a random distribution, 
Eq. | |2U . For the azimuth angle if we always used a random flat dis- 
tribution, p{(f) — const.. For the mass, we always used a flat prior 
p(M2oo) = const., whereas the a priori PDF for the concentration 
was flat in the range < C200 ^ 30 and null otherwise. The c(A/) 
relation, see Eq. ( |13t , was either enforced or neglected. Posterior 
PDFs were computed by running four chains for each of the three 
cases (only WL, only SL, combined WL and SL). We computed at 
least 20000 samples per chain and eventually added ten thousands 
more until the convergence criterium was satisfied. 

5.1 Weak lensing 

One of the main sources of bias in the lensing-based c(M) mea- 
surements is that an average triaxial halo with the major axis ori- 
ented along the line of sight appears to be over-con centrated if 
analysed under the assumpti on of spherical symmetry dOguri et al.l 
I2OO5I : ICorless & King||2007 h. Allowing for triaxial shapes and ca- 
sual orientation, the inferred distribution of either mass or concen- 
tration will come out much shallower than unde r the sharp hypoth- 
esis of spherical symmetry dCorless et al.ll2009l) . 

^ We are neglecting the third angle ip that fixes the direction in the plane 
of the sky and is related to the orientation angle 8^ 
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Tabl e 2. Intrinsic param eters inferred thi'ough lensing inversion assuming different piiors. fi and cr are the biweight estimators for location and scale of the 
PDF iBeers et aljl990l) : 71 is the skewness. Masses are in units of 10^*Mq. 



We retrieve these trends. Results are summarised in Table |2] 
In Fig. [3] we plot as binned histograms the marginalised posterior 
probability density functions for the intrinsic halo parameters, com- 
puted exploiting different a priori hypotheses. For each parameter, 
we compared the inferred PDF with the theoretical prediction from 
A'^-body simulations, see Sec. [3] Predictions for concentration and 
axial ratios were computed by convolution of A^-body predictions 
in Eqs. ( I12II13I141 andll6t with the inferred posterior PDF for the 
mass, plotted in the leftmost row. The expected and inferred C200 
distributions have different meanings. The theoretical one is the 
concentration function of a population of clusters with a mass func- 
tion as in the leftmost panel; the inferred one refers to the proba- 
bility density for the concentration of a single cluster (A1689) with 
a given mass (to be determined). The biased distribution for the 
orientation angle 9 was plotted as a continuos thick line (the same 
representation used for the random distribution of 0). The flat dis- 
tributions for qi and §2 and the random one for 6 were plotted as 
thin lines for comparison. 

The assumptions of either mildly triaxial halos, as suggested 
by simulations, or biased orientations have similar effects on the 
concentrations. The tail of the distribution at large values gets thin- 
ner whereas the maximum shifts towards lower values. Only by 
enforcing the c(A/) relation, the C200 distribution becomes very 
similar to the theoretical one, even if over-concentrated values are 
still slightly preferred. 

Since WL is not effective in determining the projected ellip- 
ticity, constraints on the intrinsic axial ratios are not sharp and are 
dominated by the a priori hypothesis employed. However, a very 
elongated shape, qi < 0.2, is strongly disfavoured with respect to 
mildly triaxial configurations with qi > 0.3. Assuming flat distri- 
butions, the mode of the distribution for gi is in the range 0.7 - 
0.8, whereas using priors from A^-body simulations, the inferred 



distributions become very similar to the theoretical ones, with peak 
values qi ~ 0.4 — 0.5. 

The results on the second axial ratio 52 and the angle (j) are 
not very significant since they are dominated by the prior. Results 
on the inclination angle of the major axis, 9, are more interest- 
ing. When assuming random orientations, aligned configurations 
(cos 6 < 1) are not significantly preferred. The peak of the PDF 
moves towards smaller values of 9 if either mildly triaxial shapes 
are enforced or the c(A/) relation is assumed to be valid. 



5.2 Strong lensing 

Strong lensing inversion was performed by deprojecting the central 
NFW halo found in the analysis of the multiple image systems. The 
likelihood was computed as £sl oc exp {— Xsl/2}. Marginalised 
posterior probability density functions for intrinsic parameters are 
plotted in Fig.|4j corresponding momenta are listed in Table [2] As 
far as the effects of priors are concerned, we retrieve the same 
trends found in the weak lensing analysis. The cluster appears to 
be over-concentrated but still compatible with predictions. Assum- 
ing a priori either mildly triaxial axial ratios or biased orientations, 
the disagreement gets smaller. 

Strong lensing gives the more effective constraint on the pro- 
jected ellipticity. Triaxiality increases the lensing strength, so halo 
models inferred from SL do not have to be extremely massive 
or concentrated. PDFs for mass and concentration obtained with 
strong lensing inversion are compatible but more peaked and less 
shallow than distributions from weak lensing. The concentration 
PDFs peak at lesser values. Central momenta for C200 shift from 
the range 13 — 6 to values between 6 — 5; PDF dispersions ob- 
tained with SL inversion are nearly one half of those in the WL 
case (1 — 2 compared to 2 — 4). There are still long tails at large 
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Figure 3. PDFs for the intrinsic halo parameters as derived from weak lensing inversion. Panels from the left to the right are for M200 , C200 , gi , <?2 , cos 9 and 
cos ip, respectively. The white binned PDF is the posterior PDF; the dark grey histogram is the distribution expected according to predictions from A'^-body 
simulations given the mass distribution in the left panel; intermediate grey bins are shared by the two plotted distributions. Thin Unes in the q'-panels denote a 
flat g-distribution. Thin lines in the 5 and (ji-angle panel denote a random orientation. The thick Une in the cos S-panel represents a biased orientation. Priors 
are stated above each set of PDFs. From the top to the bottom row, we plot the weak lensing inversion results vmder the prior hypotheses of: flat g-distribution 
and random orientation angles; flat g-distribution and biased orientation angles; flat g-distribution, random orientation angles and fulfilled c(M) relation; 
A''-body g-distribution and random orientation angles; W-body g-distribution and biased orientation angles; W-body g-distribution, biased orientation angles 
and enforced c(M) relation, respectively. Masses are in units of 10^'' Mq. 



values, but they do not extend behind C200 ^ 15. Strong lensing 
also prefers larger mass values, AM200 ~ (4 — 6) x 10^* Mq for 
the central momenta with respect to WL. The shift in mass is even 
more pronounced when the prior on the c(M) relation is enforced. 

The prior on the c(A/) relation affects the weak or the strong 
lensing inversion to a different degree. Since weak lensing PDFs 
span large ranges in both mass and concentration, the c(M) prior 
can dominate the posterior result by selecting a region which is al- 
ways compatible with the likelihood. Then, the shapes of the PDFs 
for mass and concentration with or without the prior on the (:(M) 
relation are very different. On the other hand, the concentration dis- 



tribution from SL is already peaked and in good agreement with 

theoretical predictions. Then, any c(M) prior alleviates the over- 
concentration problem for the strong lensing analysis, favouring 
lesser values of C200 and smaller dispersions, but does not dramati- 
cally change the PDFs. 

Differently from any shear analysis in the outskirts, strong 
lensing can precisely estimate the projected ellipticity in the in- 
ner regions. Whatever the priors on the axial ratios distributions, 
spherical shapes are quite disregarded. Even assuming a priori flat 
distributions, values of gi > 0.8 are now excluded and final dis- 
tributions are quite peaked favouring triaxial shapes. Using strong 
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Figure 4. PDFs for the intrinsic halo parameters as derived from strong lensing inversion. Notation and style follow Fig.|3] 



instead of weak lensing, the central momenta of the PDF for qi 
and q2 shift by ~ 0.2 and ~ 0.1, respectively. As before, when we 
adopt A'^-body like prescriptions, posterior PDFs mimic the prior 



live constraint on the projected ellipticity, the posterior distributions 
for the axial ratios are quite similar to the strong lensing results. 



5.3 Weak plus strong lensing 

In order to perform a weak+strong lensing inversion, we considered 
a combined likelihood Cam oc /Iwl x £sl. Since posterior distri- 
butions from either strong or weak lensing analyses are compatible, 
the combined investigation provides a meaningful compromise. 

Concentration values are a bit larger but still compatible with 
theoretical predictions. With respect to the inversion using only 
strong lensing, lesser values of the mass are preferred, — AA/200 ~ 
(4 — 6) X 10^ Mq, with smaller dispersions. On the other hand, 
concentrations are slightly higher, Ac2oo ~ 1 — 1.5, but with simi- 
lar dispersions and shapes. Since strong lensing gives a more effec- 



6 RESULTS 

Several features of the infeiTed distributions of intrinsic parame- 
ters are more easily understood by considering multi-dimensional 
posterior PDFs. Some trends are quite clear 



6.1 c(M) relation 

Whatever the assumptions on either orientation or shape, there 
is agreement between theoretical predictions and measured val- 
ues of mass and concentration. This can be seen analysing the 
marginalised p(M2oo, C200), see Figs.|6]and|7] where we compared 



© 0000 RAS, MNRAS 000. 000-000 



10 Sereno & Umetsu 



No c[M) prior 



Flat g-diEtribiition 



RaDdom oiieiitattons 




No r.[A'!) prior 



Flat (j-distrihiition 



mTrrmrm'i 




10 li ?D 15 




mi 




OC 02 04 OS OS IS 00 62 0,4 0.6 Oi 10 



Riased orienfatios 

4 



04 Oi 03 iD 00 0.2 04 0.6 0^ 10 



EiiToiccd c{M) 



Fiat y-diislribiiLioii 



Rimdum oiicutaliuu 



Jrrn n hrrrFfnTll 



00 02 04 0^ 03 la 00 02 04 0.6 O.F 10 



No c(M) prior 



iV -body q-distribntion 



Random oiicntattorts 



m 11 




/ 



/ 



1 



OC 02 04 Oi 03 13 00 0.2 04 0.6 OS 1.0 



No c{AI) prior 



A'^ body q distribution 



Biased orientations 



S ifl IS K) 



J I (. i B a H ■ «.2 n-l IM. "tii l.t 




i.i op 03 01 f/e OS 13 M (12 dl 06 OB 10 



Eniorced c.{M ] 



iV -body g-distribution 



Biased orientations 



S 10 :2 M iS 1! » H 



i e to iz 





1 




t — 



o lri-hi-r-r-rrr-r1 




*a 61 0.4 06 OS 13 00 0.2 04 0.6 0.3 10 



Figure 5. PDFs for the intrinsic halo parameters as derived from combined weak and strong lensing inversion. Notation and style follow Fig.|3] 



the computed posterior PDF with the theoretical conditional prob- 
ability for the concentration given a halo mass. 

A 1689 appears to be a quite typical massive cluster with a 
concentration in agreement with the tail at large values of the 
expected population of clusters of that given mass. The peak of 
p(M2oo, C200) as inferred from strong (weak) lensing inversion is 
< 1 (2)(T away from the maximum of the cluster concentration 
function, with the la tail at smaller values compatible with the me- 
dian expected value. This results is independent on the assumed 
prior. In fact, Figs.|6]and|7]consider two quite different scenarios: 
a priori flat distributions and random orientation versus axial ratios 
and inclination predicted by TV-body simulations. 

The combined weak plus strong lensing analysis further con- 
strains the preferable region in the A/200 -C200 plane. The tail at 
low concentrations is thinned, but there is still agreement. Indepen- 
dently of the priors, the inferred C200 are only > la away from the 
predicted median value. Priors from A'^-body simulations also help 



to put an upper bound on the concentrations. The 3a region for the 
inferred j3(M2oo, C200) in the weak plus strong lensing is within 

C200 < 11. 



As theoretical prediction w e considered result s from the anal- 
ysis of a full sample of clusters jPuffv et"^l2008l) . The c(M) re- 
lation for a sample of only relaxed halos has a peak at larger con- 
centrations and a smaller dispersion at a given mass. The inclusion 
of unrelaxed haloes adds a tail towards low concentrations. Even 
comparing the A1689 concentration and mass with results for re- 
laxed clusters, we would have found agreement as well. Two oppo- 
site effects counterbalance if we consider the relaxed instead of the 
full sample. On one side, the distance between the peaks of either 
the inferred or the expected distribution gets smaller; on the other 
hand, the predicted tail al large values is smaller too, reducing the 
distribution overlapping. 
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Weak lensing Strong lensing Weak and Strong lensing 




10 20 30 40 10 20 30 40 50 10 15 20 25 30 

M200 [10"M^] M200 [iO"M„] M200 [lO'-'M^] 



Figure 6. Contour plot of the marginalised PDF for Af2oo and C200 as derived under the prior assumptions of flat g-distribution and random orientation 
angles. Panels from left to right are for the weak, strong and combined lensing inversion, respectively. Contours ai'e plotted at fraction values cxp(— 2.3/2), 
exp(— 6.17/2), and exp(— 11.8/2) of the maximum, which denote confidence limit regions of 1, 2 and 3cr in a maximum likelihood investigation, respec- 
tively. The full, long-dashed and dashed lines enclose the 1, 2 and Scr regions for the predicted conditional probability c{M), respectively. 



Weak lensing Strong lensing Weak and Strong lensing 




5 10 15 20 25 30 35 10 15 20 25 30 6 8 10 12 14 16 18 20 

Mzm [W"M„] M200 [lO'^M^] M200 [lO^M^] 



Figure 7. The same as Fig. |6] when assuming an A^-body like g-distribution and biased orientation angles. 
Weak lensing Strong lensing Weak and Strong lensing 




0.0 0.2 0.4 0.6 as 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

11 91 91 

Figure 8. Contour plot of the marginalised PDF for qi and C200 as derived under the prior assumptions of flat q-distribution and random orientation angles. 
Shadowed contours and panels aie as in Fig.|6] 
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Weak lensing 




Strong lensing 




Weak and Strong lensin] 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

qi gi qi 

Figure 9. Contour plot of the marginalised PDF p{qi,q2) as derived under the prior assumptions of flat g-distribution and random orientation angles. 
Shadowed contours and panels are as in Fig.|6] The thick full, long dashed and dashed lines hmit the 1, 2 and Scr confidence regions for the A^-body like 
expected distributions. 



6.2 Shape 

Even if lensing constraints are not enough to precisely determine 
the intrinsic axial ratios, mildly triaxial halos with gi ~ 0.5 make 
a better job in fitting data than nearly spherical lenses. To see 
what can be inferred on the shape, let us focus on the results ob- 
tained using the less informative priors, i.e. flat distributions for 
the axial ratios and random orientations. In Fig. |8] we plotted 
the marginalised posterior PDF p(q\ , C200) under such hypotheses. 
Values of 0.4 < gi < 0.8 are more likely than either extremely tri- 
axial or nearly spherical shapes. Results from either weak or strong 
lensing inversion show correlation between shape and concentra- 
tion. For intermediate qi values, the dispersion in C200 is quite 
small. 

Triaxial shapes predicted by A'^-body simulations are in good 
agreement with our results. In Fig.|9] we compared shapes inferred 
from data with theoretical predictions. We used the flat prior. Even 
if a priori each configuration qi is equally likely, the posterior PDF 
prefers mildly triaxial axial ratios. Due to the better determination 
of the projected ellipticity, constraints from strong lensing are more 
significant. Axial ratios inferred in the weak plus strong lensing 
case are compatible at la confidence level with predictions from 
A'^-body simulations. Furthermore, we can exclude nearly spherical 
shapes (qi ~ g2 ~ 1) at the 3a confidence level. 



6.3 Orientation 

There are some indications for an orientations bias. Smaller values 
of 6 are preferred, pointing out to the cluster being somewhat elon- 
gated along the line of sight. This can be seen in Fig.[TOl where the 
marginalised p{cos6, C200) is plotted under the assumptions of an 
AT-body like g-distribution but for random orientations. Whereas 
the tail is extended and small values of cos 6 are still compatible, 
the peak of the distribution is at large values of cos 9. The larger the 
elongation bias, the smaller the values of the concentration compat- 
ible with data. 

We retrieved such trend for either weak (left panel in in 
Fig. fTOt and strong lensing (central panel) inversion. By combin- 
ing weak to strong lensing data, the preferred values of the con- 
centration shift towards larger values than in the only SL case, but 
the preference for small inclinations of the major axis, i.e., small 
values of 9, remains unaffected. 





Priors 










Qi 


angles 


c(M) 


WL 


SL 


WL+SL 


flat 


random 


X 


l.litO.3 


1.5±0.7 


1.5±0.7 


flat 


bias 


X 


0.9 ±0.1 


0.9 ±0.2 


1.0±0.2 


flat 


random 


/ 


1.1±0.2 


1.1±0.2 


1.0±0.3 


Af-body 


random 


X 


1.3±0.5 


1.2±0.4 


1.5±0.4 


iV-body 


bias 


X 


0.8±0.2 


0.8 ±0.2 


0.9 ±0.2 


AT-body 


bias 


/ 


0.7±0.2 


0.7±0.2 


0.8 ±0.2 



Table 3. Elongation parameter inferred by assuming different priors. Bi- 
weight estimators are reported. 



6.4 Elongation 

The cluster length along the lin e of sight can be parametrized by 
the elongation ba l lSerenoll2007h . With no regard to either the prior 
assumptions or the data-set employed (either strong or weak lens- 
ing), results are quite consistent, see Table|3] Very elongated lenses 
(ga ^ 1) may be quite strong lenses. Then, when we enforce a 
biased orientation, clusters turn out to be less massive and concen- 
trated since lensing strength is provided by the elongation. 

For halos with an orientation bias and intrinsic ratios as from 
A'^-body simulations, the line of sight elongation turns out to be 
~ 20 — 30 per cent larger than the projected dimension in the plane 
of the sky. However, due to the large statistical uncertainty, more 
spherical configurations are compatible too. When using triaxial 
axis ratios as inferred from TV-body simulations, an unrealistic tail 
at very large values for the PDF of the elongation is cut off. In the 
case of a flat prior for the qi, this tail pushes the central estimate for 
the elongation to higher values. 

These new results from lensing for the matter halo agree with a 
previous direct measurement for the gas distribut ion (eA = 1.1 ± 
0.2) obtained by combin ing X-ray and SZ data jDe Filippis et al.l 
l2005l : ISereno et al.l2006l) . The matter distribution is expected to be 
more triaxial than the intra-cluster medium, but the large errors do 
not allow us to confirm this predictio n. 

Recently, iMorandi et al.l ( 1201 ih combined lensing and X-ray 
data to infer a quite elongated halo. Assuming the cluster to be 
aligned with the line of sight and a peculiar model l ing fo r the devia- 
tion from hydrostatic equilibrium, iMorandi et al.l (1201 ll) calculated 
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eA ~ 0.6. This estimate strongly relies on the hypothesis on the 
hydrostatic status and on results from numerical simulations relat- 
ing the shapes of gravitational potential and the matter distribution, 
whic h are affected by quite large errors for very elongated clusters 
JEee & Suta2003.) . 



7 CONCLUSIONS 

Most of the tension on A 1689 as an extremely over-concentrated 
lensing cluster can be eliminated. We exploited two main ingre- 
dients: i) a full triaxial investigation of the lensing halo; ii) a 
Bayesian statistical analysis. Simulations and theoretical analyses 
show that massive halos are quite triaxial. A rich variety of shapes 
and orientations has to be considered, some of which allow moder- 
ately concentrated halos to be lenses as effective as much more in- 
trinsically concentrated clusters. In particular, lensing cross-section 
for ellipsoidal halos elongated along the line of sight are enhanced. 

Bayesian methods allow us to explore the full parameter 
space, avoiding to restrict the analysis in the vicinity of the con- 
figuration with the maximum likelihood. A proper investigation of 
the parameter space provides also a more reliable estimate for the 
statistical uncertainty. 

Our results are in good agreement with predictions based on 
numerical simulations. The triaxial lensing analysis of A 1689 sug- 
gests a mildly triaxial cluster. Computed concentrations (C200 ~ 
10 ± 3 or 7 ± 1 for the only WL or the combined WL-l-SL case) are 
only slightly larger (< 2a) than the median concentration (C200 — 
3 . 2 ± 1 . 6) expected for the cluster mass M200 = 13 xlO"M0. In- 
versions of either strong or weak lensing data are reassuringly com- 
patible, and the combined analysis enhanced the picture of A1689 
as a typical massive lensing cluster. There is some evidence for an 
alignment of the major axis with the line of sight. The size of the 
matter distribution along the line of sight is larger by ~ 20 per cent 
than the width in the plane of the sky, but more spherical configu- 
rations are still compatible with data. 

T he theoretical form of the c(M) relation is still debated. Re- 
cently, IPrada et aEI ( 1201 ll) showed that the dependence of concen- 
tration on halo mass and its evolution can be obtained from the 
root-mean-square (rms) fluctuation amplitude of the linear density 
field ij{M, z). They noticed a flattening and upturn of the relation 
with increasing mass. Estimated concentrations for galaxy clusters 
might be substantially larger than previous results. Agreement with 



our result for A 1689 would consequently be even better, since they 
predict a median C200 ~ 5.1 for a cluster of A/200 = 13xlO'^*M0 
at the redshift of A1689. Accurate determinations of n^, the spec- 
tral index of the primordial power spectrum, and as, the rms am- 
plitude of linear mass fluctuations, are also needed for reliable the- 
oretical prediction on c(M). 

Due to the different statistical approaches, it is difficult to 
compare our estimates with previous works. Our analysis is in 
broad agreement with previous results which found a concentra- 
tion of C200 ~ 6 for strong lensing analyses or ~ 8-12 for weak 
lensing ana lyses. Previous result s from literat ure have al ready been 
collected in fcorless efZI ( l2009l table 4) and lCoe et ^ Jloia ta- 
ble 2). In general, triaxiality allows for smaller conc entrations and 
larger dispersions. Assuming spherical symmetry, lUmetsu et al.l 
('2009) found C200 ~ 12. 5l^:? from their one-dimensional WL 
analysis of Subaru distortion data. Due to the strict a priori hy- 
pothesis on the intrinsic shape, this estimate is higher than what 
obtained in the present paper allowing for triaxiality. The disper- 
sion for the spherical symmetric case is smaller than our WL result 
assuming a flat prior on axial ratios and random orientations. Forc- 
ing the shape to follow results from A'^-body simulations and the 
orientation to be biased, dispersion is comparable with the spheri- 
cal symmetric case, but the estimated concentration is significantly 
lower (C200 ~ 10 vs. C200 ~ 13). 

The most m eaningful comparison can be made with 
ICorless et ai] ( l2009l) . who first employed a Bayesian analysis to de- 
project weak lensing data of A 168^0 Even if they used a different 
data set obtained with the ESO/MPG Wide Field Imager, the ob- 
tained posterior distributions for C200 are in agreement with ours 
for both location and dispersion. 

Orientation and shape biases can be very sizable when esti- 
mating halo concentrations from lensing but are less effective in 
X-ray analyses (Oguri & Blandford 2009; Meneghetti e t al. 201(t) 
which usually provide lesser values for C200 in A 1689 ICoe et al.l 
I2OIOI) . A proper statistical modelling of triaxiality is then manda- 
tory for reliable estimates. This deserves attention since an unbi- 
ased lensing estimate can provide concentration values without as- 
suming hydrostatic equilibrium, which can plague X-ray methods 



^ The "flat" model in ICorlessetai]b009h used priors similar to our flat q- 
distribution with random orientations; their "Shaw" model used priors for 
the axial ratios from numerical simulations. 
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( iMolnar et allllOlOl) . Combined X-ray and SZ methods are inter- 
esting too, since thiey can directly infer thie elongation of tlie intra- 
cluster m edium distribution without the hydrostatic equilibrium hy - 
pothesis jPe Filippis et al.ll2005l : [Sereno et al.l200^ : ISerenoll2007h . 
However, the gas density is expected to differ from the dark matter 
profile. Some further hypotheses have to be used to link the gas to 
the dark matter distribution. 
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